Quasiparticles governing the zero-temperature dynamics of the ID spin-1/2 
Heisenberg antiferromagnet in a magnetic field 
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The T — dynamical properties of the one-dimensional (ID) s — | Heisenberg antiferromagnet 
in a uniform magnetic field are studied via Bethe ansatz for cyclic chains of N sites. The ground 
state at magnetization < M z < N/2, which can be interpreted as a state with 2M Z spinons or as 
a state of M z magnons, is reconfigured here as the vacuum for a different species of quasiparticles, 
the psinons and antipsinons. We investigate three kinds of quantum fluctuations, namely the spin 
fluctuations parallel and perpendicular to the direction of the applied magnetic field and the dimer 
fluctuations. The dynamically dominant excitation spectra are found to be sets of collective excita- 
tions composed of two quasiparticles excited from the psinon vacuum in different configurations. The 
Bethe ansatz provides a framework for (i) the characterization of the new quasiparticles in relation 
to the more familiar spinons and magnons, (ii) the calculation of spectral boundaries and densities of 
states for each continuum, (iii) the calculation of transition rates between the ground state and the 
dynamically dominant collective excitations, (iv) the prediction of lineshapes for dynamic structure 
factors relevant for experiments performed on a variety of quasi-lD antiferromagnetic compounds, 
including KCuF 3 , Cu(C 4 H 4 N2)(N03)2, and CuGe0 3 . 



I. INTRODUCTION 

Quantum spin chains are some of the most intensively 
studied models representing strongly fluctuating quan- 
tum many-body systems because of their amenability 
to exact analysis and because of the sustained interest 
in materials exhibiting quasi-one-dimensional magnetic 
properties. Of particular interest are the dynamical prop- 
erties in the low-temperature regime, reflecting strong 
quantum fluctuations. 

Quantum fluctuations result from the time evolution of 
nonstationary observables of a many-body system in the 
ground state. They can be investigated (experimentally, 
theoretically, or computationally) by dynamical probes. 
The three main ingredients of each dynamical probe, (i) 
interaction Hamiltonian, (ii) ground state, (iii) dynami- 
cal variable, make a specific set of collective excitations 
visible to the probe. The specificity is determined by the 
symmetries of all three ingredients. 

A dynamical probe yields information on spectrum and 
transition rates. Different sets of data are collected from 
the same many-body system [ingredients (i) and (ii)] via 
particular fluctuation operators [ingredient (iii)]. Differ- 
ent views of the quantum fluctuations are filtered out, by 
operator specific selection rules and transition rates.EJ 

Collective excitations are modes in which some of the 
tightly coupled fundamental degrees of freedom (elec- 
trons, ions, atoms) move collectively in more or less 
complex patterns. The free-particle like normal modes 
known to exist in systems made of linearly coupled de- 
grees of freedom are the inspiration of attempts to de- 
scribe collective excitations quite generally as composites 
of elementary modes that are weakly coupled at most. 
This requires that the ground state of the system can be 
meaningfully interpreted as a physical vacuum in which 
certain kinds of elementary excitations (quasiparticles) 



move without attenuation and scatter off each other non- 
destructively. 

In completely integrable many-body systems, the iden- 
tity of the quasiparticles in any given cigenstate is upheld 
on a rigorous basis and encoded by a set of quantum num- 
bers. All excited states can then indeed be generated 
systematically via the creation of quasiparticles from the 
ground state configured as a physical vacuum. The in- 
teraction of the quasiparticles may not be weak, but it is 
of a kind which preserves their identity. The factorizabil- 
ity of the associated S'-matrices, which is characteristic 
of completely integrable systems, reduces all quasiparti- 
cle couplings to two-body scattering events for which a 
general solution can be. formulated, e.g. in the form of a 
Bethe wave function BlI 

The focus here is set on the quasiparticles which gov- 
ern the quantum fluctuations of the one-dimensional (ID) 
s = | Heisenberg antiferromagnet in an external mag- 
netic fieldfli 



A" 



H = [ JS » • S «+l - hS n] 



(1) 



The ground state at h > h s = 2 J, \F) = | ft • • • T), has 
saturated magnetization, M z = N/2. It is the reference 
state of the coordinate Bethe ansatz and plays the role 
of the vacuum for magnons (spin-1 quasiparticles). All 
eigenstates of H are described as configurations of inter- 
acting magnons. The ground state at h = 0, \A), has 
magnetization M z = 0. It contains N/2 magnons. In the 
framework of the Bethe ansatz, it is reconfigured as the 
physical vacuum for spinons, a species of spin-^ quasi- 
particles, and the entire spectrum of H is reinterpreted 
as composites of interacting spinon pairs. Likewise at 
intermediate fields, < h < hs, the ground state \G) 
is reconfigured as a new physical vacuum, and the low- 
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lying collective excitations are most effectively described 
as composites of two new species of quasiparticles, named 
psinon and antipsinon. 

In a recent paper a detailed description of these 
quasiparticles in the framework of the coordinate Bethe 
ansatz was given. Their role in the zero-temperature spin 
fluctuations parallel to the direction of the magnetic field 
was elucidated in the form of lineshape predictions for the 
associated dynamic structure factor. Here we present a 
more comprehensive set of applications, which also in- 
cludes the perpendicular spin fluctuations and the dimer 
fluctuations. 

Physical realizations of Hciscnberg antifcrromagnetic 
chains have been known for many years in the form of 3D 
crystalline compounds with quasi-lD exchange coupling 
between magnetic ions. For the study of magnetic-field 
effects in the dynamics as predicted in this paper, the 
coupling must not be too weak or else it will be hard to 
reach the low-temperature regime. It must not be too 
strong either or else it will be hard to reach a magnetic 
field that makes the Zeeman energy comparable to the 
exchange energy. One compound that promises to be 
particularly suitable for this purpose is copper pyrazine 
dinitrate [Cu(C 4 H 4 N 2 )(N03)2]-l3 

The spin fluctuations can be observed directly via in- 
clastic neutron scattering experiments. At very low tem- 
peratures, the dominant transitions in the scattering ex- 
periment are between the ground state |G) and a set of 
excitations |A) that are reachable by one of the spin fluc- 
tuation operators Sjf — N~ 1 / 2 J2 n e%qn S£, [i = x,y,z. 
In the T = dynamic spin structure factors 

S^( q ,Lo) = 2nY / \(G\SZ\\)\ 2 6(cu~Lu x ), (2) 

A 

each transition with lu\ = E\ — Eq and q = k\ — fcg 
contributes a spectral line of intensity 2ir\ (G\S£\\) \ 2 . 

Some quasi-lD antiferromagnetic compounds, of which 
CuGe03 is the most prominent example |ml!iJ are sus- 
ceptible to a spin-Peierls transition, which involves a lat- 
tice distortion accompanied by an exchange dimerization. 
The dimer fluctuations, D q = N~ 1/2 £ n e lqn S n ■ S n+1 , 
as captured by the dynamic dimer structure factor 

S DD (q,u) = 2nY / \(G\D q \\)\ 2 5 - u x ) (3) 

A 

may not be as directly observable as the spin fluctuations 
but an understanding of their quasiparticle composition 
is a matter of no less importance. 

II. MAGNONS, SPINONS, PSINONS 

The coordinate Bethe ansatz provides a natural clas- 
sification of the eigenstates of ([!]) in terms of interacting 
magnons. The structure of the Bethe wave function, its 
determination via the solution of the Bethe ansatz equa- 
tions, and its use for the calculation of matrix elements 
are summarized in the Appendix. 



For our discussion here it turns out to be sufficient to 
consider r-magnon scattering states of the set K r . In 
the invariant Hilbert subspace of magnetization M z — 
N/2—r, the Bethe quantum numbers of this set comprise, 
for < r < N/2 and < to < N/2 — r, all configurations 

~ + ^-m<h<h<---<I r <^- -+m. (4) 

The Bethe ansatz suggests a threefold interpretation of 
the ground state \G) at < M z < N/2 with quantum 
numbers 

r n J N M 1 N M z 1\ 

Depending on the reference state (pseudo-vacuum) used, 
it can be regarded as a scattering state of N/2 — M z 
magnons, a scattering state of 2M Z spinons, or the phys- 
ical vacuum of psinons II 

The states in the set K r then all contain the same 
number of magnons or spinons but different numbers 
of psinons. The integer quantum number to selects all 
states from K r that contain to pairs of psinons. The 
ground state \G) at M z = N/2 — r is the only state with 
to = 0. The quasiparticle role of the psinons in the 2- 
psinon (to = 1) and 4-psirupn (to = 2) scattering states 
was highlighted previously^ 

The excitations which are important in S zz (q, u>) (par- 
allel spin fluctuations) at M z = N/2 — r were found 
to consist of a small subset of K r which includes 2to- 
psinon states over the entire range to. However, all 2to- 
psinon states with significant spectral weight were found 
to belong to particular configurations of Bethe quantum 
numbers 7, in which 2m — 1 psinons behave like a sin- 
gle degree of freedom with properties akin to those at- 
tributed to an antiparticle. The spectrum of S zz (q,u>) 
was thus identified as arising predominantly from psinon- 
antipsinon (ipip*) excitations Jj Here our goal is to iden- 
tify and interpret the dynamically relevant excitations 
also for S xx (q,oj) = j[S + -(q,uj) + S- + (q,uj)} (perpen- 
dicular spin fluctuations) and Sdd(q,w) (dimer fluctua- 
tions), where we expect psinons and antipsinons to occur 
in different combinations. 



III. SYMMETRIES AND CONSEQUENCES 

Narrowing down the dynamically dominant sets of ex- 
citations and characterizing them as specific quasiparti- 
cle configurations proceeds in three steps. First we limit 
the set of relevant excitations by the application of selec- 
tion rules which are imposed by the symmetry proper- 
ties of the Hamiltonian ([l]) and the fluctuation operators 
Sg, Sq, D q and which are valid for arbitrary system sizes. 
Then we identify additional selection rules that are valid 
only for N — > oo. Finally, we identify from the states not 
yet excluded those whose transition rates are predomi- 
nant in S zz (q,uj), S-+(q,u), S+-(q,w), and S DD (q,uj). 
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This last step, which here is carried out empirically, may 
very well find its ultimate justification bjs further sym- 
metries related to complete integrabilityJlil 



A. Selection rules for arbitrary TV 



The conservation laws of the total spin St and its z- 
component S T imply that transitions between eigenstates 
of (|l|) induced by the (nonstationary) spin fluctuation op- 
erators S*, Sg (vector) and the dimer fluctuation oper- 
ator D q (scalar) satisfy stringent selection rules. The six 
classes of excitations with permissible transitions from 
\G) with St — S T — M z for the fluctuation operators 
Sg , S^ , D q are all listed in Table |[ The locations of these 
classes of excitations relative to the ground state in the 
(St, S^)-plane are shown in Fig. |l]. 
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FIG. 1: Transitions between the ground state \G) with quan- 
tum numbers St = S^ = M z and six classes of excitations 
permitted by microscopic selection rules (Ref. [Hi]) . Each class 
may contribute to exactly one of three dynamic spin struc- 
ture factors. Class (ii) also contributes to the dynamic dimer 
structure factor. Fewer classes of permissible excitations exist 
for M z = (h = 0) and M z = TV/2 (h = h s ). 



B. Selection rules for TV — > oo 

Before we begin evaluating matrix elements from Bcthc 
ansatz solutions in production mode, we take note that 
the rotational symmetry of the Hamiltonian (Q) and the 
vector nature of the spin fluctuation operator (Sg , 5|, Sg) 
imply the following rigorous relations between transition 
rates invoicing excitations that belong to the same St~ 
multipletJl3 



\(G\S Z g\\V)f 

\(G\S+\\W)f 
\{G\S+\\M)f 



\(G\s-\\m)\* 

2(M Z + 1) ' 
|(G|5-|A('»))| 2 

(M z + l)(2M z + \y 

2|(G|^|A(»>)| 2 
M z 



(6a) 
(6b) 
(6c) 



The significance of the relations (||) is not limited to 
their usefulness in reducing computational work. The 
magnetization is an extensive quantity, implying M z cx TV 
at h 7^ 0. All transition rates for class (i), (iv), (v) exci- 
tations are then suppressed by factors TV or TV 2 relative 
to the transition rates of class (ii), (iii) excitations. The 
consequence is that in a macroscopic system at h ^ 0, 
the spectral weight of all class (i), (iv), (v) excitations 
in the T =p-0 dynamic spin structure factors S flfl (q,uj) is 
negligible .E3 



C. Selection rules related to integrability 

We shall find empirically that in each one of the re- 
maining classes (ii), (hi), and (vi) there exists a two- 
parameter set of excitations which governs one of the dy- 
namic structure factors of interest here. We shall name 
these sets P2, P3, and P6, respectively (see Table |). 
Corresponding finite- A spectral contributions of class (i) , 
(iv), and (v) states can be inferred from Eqs. (^|). The 
associated sets PI, P4, and P5 have their position in the 
(q, w)-plane shifted vertically relative to the sets P2, P3, 
and P6 because of the Zeeman splitting, and the spectral 
weight of the former is suppressed by factors A or A 2 as 
explained previously. 



IV. DYNAMICALLY DOMINANT 
EXCITATIONS 



We note that classes (ii), (iii), and (vi) include the 
sets K r for r = R, R — 1, R + 1, respectively, while the 
remaining three classes (i), (iv), and (v) include sets of 
states that belong to the same St multiplets as the sets 
Kr or Kr+x. Table | also lists the Bethe quantum num- 
bers = 1, . . . , r for the type-A r or equivalent states 
and describes how the rapidities z,, i = 1, . . . , r for these 
states are obtained. 



In our search for the dynamically most relevant excita- 
tions, we focus on the case of magnetization M z /N = \ 
(half the saturation value). We explore the transition 
rates for the spin and dimer fluctuation operators be- 
tween the ground state \G) and the type-A r states in the 
classes (ii), (iii), (vi). These excitations are found to con- 
tribute most of the spectral weight to the dynamic spin 
and dimer structure factors. 
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TABLE I: Specifications of type-A" r or equivalent states from six classes. Each class contains states that contribute to a specific 
dynamic spin structure factor at T = 0. Class (ii) also contributes to the dynamic dimer structure factor. All specifications are 
relative to a given ground state with S T = St = N/2 — R = M z , where M z is the magnetization in a field of a certain strength 
h. The last column identifies the three subsets of excitations which dominate the spin and dimer fluctuations for N — > oo. 



class 


St 
St 


S m (q,oj) 


r 


Bethe quantum numbers 


Bethe ansatz solutions 


dynamically 
dominant sets 


(0 


M z + 1 

M z 


S zz (q,u>) 


R 


/W = / («i) + i )i = 1) ... >R _ 1 
jW = l(N-R+l) 


« i W =« i (M<) ,i = l,...,fl-l 

«jj = OO 




(ii) 


M z 
M z 


S Z z(q,u) 


R 


l\ H \i = l,... ,R 
from (p|) with r = R 


*f\i = l,... ,R 
from (A3) with r = R 


ipip* (P2) 


(iii) 


M z + 1 
M z + 1 


S- + (q,u) 


R-1 


lf"\i = l,... ,R-1 
from (§) with r = R - 1 


z? ii) ,i = l,...,R-l 
from (A3) with r = 7? - 1 


ipip (P3) 


(iv) 


M z + 1 
M z - 1 


S+~(q,u) 


R+l 


/(») =/(«*)+!, i = i,... 


zf v) =zf i \i = l,... ,R-1 

(iv) f* 1 ') 
Z R ~ Z R+± ~ 00 




(v) 


M z 

M z - 1 


S+-(q,u) 


R+l 


/w = 4^ + i,i = i,... 


z^=zf\i = l,...,R 
z rU = °o 




(vi) 


M z - 1 
M z - 1 




R+l 


J**>,i = l,... 

from (§) with r = R + 1 


z[ vi \i = l,... ,R+1 
from (A3) with r = R+l 


ipip* (P6) 



A. Perpendicular spin fluctuations (P3) 

The spectral weight in S \-(q, u>) is carried exclusively 

by class (iii) excitations (see Table |j]). A systematic 
study of the transition rates of type-iv^-i states for 
in = 0, 1, 2, . . . reveals that the dominant contributions 
to the spectral weight come from two-psinon (ipip) states. 
The Bethe quantum numbers of the states involved in 
these transitions are shown in Fig. || for N — 16 and 
M z = 4. The top row represents the ground state 
\G), which contains four magnons (small circles) or eight 
spinons (large circles). This is the psinon vacuum at 
R — 4. The next row is the lowest-lying two-psinon (tpip) 
state excited from \G). This excitation also plays the role 
of the psinon vacuum (m = 0) at R — 3. It is then char- 
acterized as containing three magnons or ten spinons. 

Mobilizing the two innermost spinons turns them into 
psinons (gray circles). The remaining five rows in Fig. ^ 
represent ipip states (m = l)at0<g<7r for R = 3. 
The state in the second row is also counted as a ipip state. 
Hence there are six of them in total for N = 16. Freeing 
up two additional spinons from the sidelines produces a 
set of four-psinon states (m = 2) of which the ipip states 
(m = 0, 1) are special members. The maximum number 
of psinons that can be mobilized at R = 3 is equal to the 
number of spinons: 2M Z = 10. 

We have calculated the transition rates |(G|5~|A^ m ))| 2 
between the ground state for N = 24, M z = 4 and all 
2m-psinon excitations. We found that the ipip states are 
predominant. They are listed in Table || along with the 
momentum and energy transfer of the associated spectral 
lines. For N — > oo the ipip states form the continuum P3 
in (q, w)-space. In Fig. ||(a) we have plotted all states 
belonging to P3 for N = 64 (circles) and the spectral 
boundaries for N — > oo. The range of P3 is restricted to 



Qs _ |<z| _ 1") where 

q s ee 2ttM z /N. (7) 

Note that the continuum P3 is displaced by Aq = tt rel- 
ative to the two-psinon continuum discussed in the con- 
text of Ref. H The psinon vacuum used for P3 is the 
state in the first row of Fig. |^. In Ref. ^ the state in the 
second row is the vacuum. These choices are dictated by 
the different fluctuation operators considered now and 
then. 

The relative integrated intensity of the ipip states, 
St^(q)/S-+(q), is plotted in Fig. |(b) for various N at 
fixed M z = N/4. Corresponding data for the absolute 

integrated intensity S ^(q) are shown in the inset. We 

observe that there is virtually no intensity at \q\ < q s , 
outside the range of continuum P3. As \q\ increases from 
q s toward tt, S increases gradually and at an ac- 
celerated rate. The value at the zone boundary diverges 

I; -6-5-4-3-2-1 1 2 3 4 5 6 q 

I I I I I I I I I I I I I 

oooo««««oooo o 

OOOOO* • •OOOOOs 
OOOOO* •0»00007 

ooooo»o» •OOOOe 

OOOOOO* • •OOOOs 
0000»0»«»00008 

oooo»oo» •oooo? 

FIG. 2: Psinon vacuum |G) for N — 16, M z = 4 and ipip states 
with q > from the set K3 out of class (iii). The Ii values 
are marked by the positions of the magnons (small circles). 
The spinons (large circles) mark /^-vacancies. A subset of the 
spinons are called psinons (gray circles) The wave numbers 
q = k — kc are given in units of 2n/N. 
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TABLE II: ipip states with q = k — ka > from the set 
Kr-\ out of class (iii) excited from the psinon vacuum \G) 
for N = 24, R — 6: Bethe quantum numbers, wave number 
(in units of 2n/N), energy, and transition rate. The ground 
state has &g = and Eg = —11.5121346862 and is realized 
at h = 1.58486 ... for N -> oo. 







2/, 






q 


E -E G + h 


\(G\S-\\ (M) )\ 2 


-2 +0 


+2 


+4 


+6 


7 


1.8184194057 


0.0626472812 


-4 +0 


+2 


+4 


+6 


8 


1.9556536638 


0.1424479327 


-4 


-2 


+2 


+4 


+6 


9 


2.0088932004 


0.2414561733 


-4 


-2 


+0 


+4 


+6 


10 


1.9746665911 


0.3895105465 


-4 


-2 


+0 


+2 


+6 


11 


1.8556790610 


0.6950450037 


-4 


-2 


+0 


+2 


+4 


12 


1.6606927007 


2.2311978057 


-6 


+0 


+2 


+4 


+6 


9 


2.1612403300 


0.0104636629 


-6 


-2 


+2 


+4 


+6 


10 


2.2113923024 


0.0348248694 


-6 


-2 


+0 


+4 


+6 


11 


2.1742818977 


0.0780182561 


-6 


-2 


+0 


+2 


+6 


12 


2.0527517142 


0.1743614335 


-6 


-4 


+2 


+4 


+6 


11 


2.3389550993 


0.0062315774 


-6 


-4 


+0 


+4 


+6 


12 


2.2987048423 


0.0260883689 



in the thermodynamic limit: S |_(7r) ~ N\ l/ 71 ., with an 

exponent rj{M z /N) that assumes the valueoEj 

77(1/4) = 1.53122... (8) 

for the situation at hand. It reflects the divergence 
S-+{q) ~ |tt - q] 1 ^- 1 for N = 00. The relative ij)ip 

contribution to S \.(q) rises rapidly from zero at q > q s 

toward a value exceeding 97.8% at q — tt. The solid 
line in the inset is obtained from a two-parameter fit, 
a\n - qf/i- 1 + 6, of the data at q > tt/2. 

When we decrease M z at fixed N, the soft mode at 
q = tt remains stationary while the soft mode at q = q s 




FIG. 3: (a) Energy versus wave number of the tpip exci- 
tations at M z = N/4 for N = 64 (circles) and AT — > 00 
(continuum P3 between curves), (b) Integrated intensity of 
S- + (q) (inset) and relative iprp contribution (main plot) for 
N = 12,16,20,24,28. 



moves to the left. At M z = 0, the ipip states become the 

two-spinon triplets. The twej-sainon part of S \-(q,w) is 

exactly known for N — > ooJl3'tZI Conversely, when we in- 
crease M z , the soft mode at q — q s moves to the right and 
thus narrows the range of P3 continually. At saturation 
(M z /N — 4) the function S y{q,u) vanishes identically. 



B. Dimer and parallel spin fluctuations (P2) 

The parallel spin fluctuations were already analyzed in 
Ref. ||. The relevant excitations are contained in the set 
Kr out of class (ii). This set is subdivided into sets of 
2m-psinon states for < m < M z = N/2 — R. Each set 
for m > contributes one branch of excitations with sig- 
nificant spectral weight to S zz {q,u). Figure ^ shows the 
configurations of Bethe quantum numbers Ii of all these 
states for N = 16, M z = 4. The top row represents the 
psinon vacuum (m = 0). The four groups of states un- 
derneath represent the dynamically dominant branches 
of 2m-psinon states for m = 1, . . . ,4. 

We argued that the ij-configurations of these excita- 
tions suggest a simpler interpretation in terms of two 
quasiparticles, namely one psinon (ip) and one antipsinon 
(if)*)- We reinterpreted the series of dynamically domi- 
nant branches taken from 2m-parameter sets of multiplc- 
psinon states as a single two-parameter set of i/ji/j* scat- 
tering states. These states form the continuum P2 for 
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FIG. 4: Psinon vacuum \G) for N = 16, M z = 4 and set of 
ipip* states with < q < n out of the set K4. The Ii are 
given by the positions of the magnons (small circles) in each 
row. The spinons (large circles) correspond to Ji-vacancies. 
The psinon (ijj) and the antipsinon (ip*) are marked by a large 
and a small grey circle, respectively. 
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TABLE III: ipip* states with q = k — ko > from the set 
Kb. out of class (ii) excited from the psinon vacuum \G) for 
TV = 16, R = 4: quantum number m, wave number (in units 
of 2-n/N), energy, and transition rates. The ground state 
has kG — and Eg — —11.5121346862 and is realized at 
h — 1.58486 . . . for TV — * oo. The Bethe quantum numbers 
were listed in Table II of Ref. ^. 
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1.2708459328 
1.4285177129 
1.3858078992 
1.1488426600 
1.6819046570 
1.8257803105 
1.7724601200 
1.5309413164 
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0.0426892576 
0.0552255878 
0.0743667351 
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0.0753564030 
0.1406415212 
0.0235815843 
0.0443726010 
0.0744641955 
0.1686893882 



1.0000000000 
0.1201967890 
0.1687346681 
0.2298023543 
0.3456324084 
0.0720507048 
0.1098585317 
0.1555227849 
0.2470269183 
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0.0866741700 
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0.0060903835 
0.0140423747 
0.0259881320 
0.0589091070 



FIG. 5: (a) Energy versus wave number the ipip* excitations 
at M z /N =1/4 for TV = 64 (circles) and TV -> oo (partially 
folded continuum P2 outlined by solid lines), (b) Integrated 
intensity S Z z(q) (inset) and relative P2 contribution (main 
plot) for TV = 12, 16, 20, 24, 28, 32. (c) Integrated intensity 
So -D (q) (inset) and relative P2 contribution (main plot) for 
TV = 12, 16, 20, 24, 28. The lines in (b) connect the TV = 32 
data points and the lines in (c) the TV = 28 data points. The 
values of S zz (tt/2) and Sdd{"k/2) extrapolated to TV — > oo 
are marked (+). 



TV — * oo in the (q, cj)-plane as illustrated in Fig. |](a). 
From the shape of the continuum with its soft modes at 
q = and q = q s , where 



7T - 2irM z /N, 



(9) 



and with the partial overlap along a stretch of the up- 
per boundary, we reconstructed the energy-momentum 
relations of the tp and ip* quasiparticlesfl 

The corresponding search for the dynamically domi- 
nant dimer excitations again points to the tpip* contin- 
uum P2. The tpip* transition rates for the fluctuation 
operators S* and D Q in a system with TV = 16, M z = 4 



arc listed in Table III for comparison. These data sug- 



gest that the spectral weight in Sdd{cl,u) is concen- 
trated more heavily at lower energy than is observed in 
S zz {q, u). A more quantitative discussion of this evidence 
will follow in Sec. |y|. 



Finite- TV data for the integrated intensities S zz (q) and 
Sdd(q) are presented in Figs. ||(b) and (c), respectively. 
Both static structure factors rise from zero at q = to a 
cusp-like maximum at q = q s = tt/2, where the soft mode 



is located. The cusp is of the form 



ii' 1 . A two- 



parameter fit, aN 1 ~ ri + b, of the data at q = q s yields the 
extrapolated values S zz (n/2) ~ 0.307 and Sdd{k/2) ~ 
0.641. On approach to q = w, the intensity drops more 
drastically in SDD(q) than in S zz (q). 

For both kinds of fluctuations, the intensity at q < q s 
is almost exclusively originating from ipip* excitations. 
At q = q s , the relative iptp* contributions to S zz (q) and 
Sdd(q) are estimated to be at least 93% and 95%, respec- 
tively, in the limit TV — > oo. At q > q s the ipip* parts of 
S zz (q) and Sdd{q) decrease monotonically but remain 
dominant except in the immediate vicinity of the zone 
boundary. The data suggest a qualitative difference in 
how the relative ipip* intensities approach zero as q — > tt. 
If the behavior near the zone boundary can be described 
by a power-law, ~ \n — g| 7 , then we predict 7 > 1 for t 
dimer fluctuations and 7 ~ 0.3 for the spin fluctuations 

Upon varying the value of M z , the continuum P2 
changes its shape continuously. In both limits, M z /N — > 
and M z /N — > |. it degenerates into a single branch 
and then vanishes .Enl3 At M z = the dimer fluctu- 
ations and the spin fluctuations are produced by en- 
tirely different sets of excitations. S zz (q,u>) is known 
to be dominated by the continuum of two-spinon triplet 
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excitations,t£l which are Bethe ansatz solutions with real 
rapidities. 5dd(<7,ijj) is presumably governed by two- 
spinon singlet excitations, which are Bethe ansatz solu- 
tions with complex rapidities J13 



C. Perpendicular spin fluctuations (P6) 

For finite N, the spectral weight in the dynamic struc- 
ture factor S-\ (g, uj) probes excitati ons fr om classes (iv)- 

(vi). However, we know from Sec. Ill B that the inten- 
sities of class- (iv) and class- (v) excitations are bound to 
fade away in the limit N — > oo. The class- (vi) excitations 
contain the set Kr+\. It turns out that much of the spec- 
tral weight in S-\ (q, u>) is carried again by tf>tp* excita- 
tions. However, the transitions (G\S^\X) probe the ipip* 
states in a different invariant subspace than the transi- 
tions (G\Sg\X), (G\D q \X) do. This causes some dramatic 
changes in the spectrum and in the spectral-weight dis- 
tribution. 

The differences are best illustrated by Fig. |^ in relation 
to Figs. || and ||. The top row in all three figures shows 
the /^-configuration of the ground state \G) for N = 16, 
M z = 4. The remaining rows in Fig. [| represent the 
ipip* states in the same invariant subspace, whereas the 
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FIG. 6: Psinon vacuum |G) for N = 16, M z = 4 and set of 
iptp* states with < q < n out of the set K5. The U are 
given by the positions of the magnons (small circles) in each 
row. The spinons (large circles) correspond to /^-vacancies. 
The psinon (ip) and the antipsinon (tp*) are marked by a large 
and a small grey circle, respectively. The last row describe a 
state from class (v) which belongs to the same SV-multiplet 
as the psinon vacuum. 
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FIG. 7: (a) Energy versus wave number the tpip* excitations 
at M z = AT/4 for N = 64 (circles) and N -» 00 (partially 
folded continuum P6 outlined by solid lines), (b) Integrated 
intensity S+-(q) (inset) and relative P6 contribution (main 
plot) for TV = 12, 16, 20, 24, 28. The solid line in the inset 
results from a two-parameter fit as explained in the context 
of Fig. § 



remaining rows in Fig. g represent the ipip states in the 
subspace with one less magnon (i.e. two more spinons). 
The first ipip state (second row in Fig. |2[) also plays the 
role of the psinon vacuum in that subspace. 

The second row in Fig. [| represents the lowest exci- 
tation probed by (G|iS+|A) and, at the same time, the 
psinon vacuum with one more magnon (i.e. two less 
spinons). The three groups of five states underneath rep- 
resent the complete set of ipip* states in the same invari- 
ant subspace. Because the momentum transfer for these 
ipip* states is relative to a different psinon vacuum than 
was the case for the ipip* states discussed in Fig. [|, the 
observable spectrum of the continuum P6 which emerges 
for N -> 00 [Fig. 0(a)] is the mirror image of the contin- 
uum P2 [Fig. |(a)]. 

The last row in Fig. |6| is not exactly a ipip* state. It 
differs from the iptp* in the previous row only by the 
smallest change in one Bethe quantum number. This 
class- (v) state belongs to the same Sr-multiplet as the 
ground state \G) (top row). We have included it here 
because its transition rate is significant. In fact, it is 
the only excitation at q = with a nonzero transition 
rate. Even though it is not a member of the set Kr + i, 

its contribution to (q, to) marks a natural endpoint of 

the continuum P6. The excitation energies and transition 
rates pertaining to all states shown in Fig. |^ are listed in 
Table [IV. 

The integrated intensity (q) as shown in Fig. (?](b) 

is almost flat in the region < q < q s . The intensity at 
q = is exactly known:E3 



S+_(0) = 2M Z /N = 1/2. 



(10) 



At q > q s the function (q) rises gradually and with 

increasing slope ending in a divergence at q = n. The 
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TABLE IV: ipip states with q = k - k G > from the set 
Kn+i out of class (vi) excited from the psinon vacuum \G) for 
N — 16, R = 4: Bethe quantum numbers, quantum number 
m, wave number (in units of 2n/N), energy, and transition 
rate. The ground state has k G = and E G = -11.5121346862 
and is realized at h — 1.58486 . . . for N — * oo. The last row 
describe a state from class (v) which belongs to the same 
Sr-multiplet as \G). 
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for very large N via Bethe ansatz. The factor M(q,ui) 
represents the scaled transition rates N\(G\S£\\)\ 2 , \i — 
z, +, — , or N\(G\D q \X)\ 2 between the ground state and 
the sets of iptjj or i/ji/j* states. These matrix elements 
are also calculated via Bethe ansatz but only for much 
smaller systems. 

For the applications considered here, the product 
ansatz depends on a reasonably fast convergence, within 
the spectral boundaries of the continua P2, P3, and P6, of 
the finite- TV transition rate data toward a smooth func- 
tion M(q, uj) as N — ► oo. Problems with this ansatz 
arise when the finite-size excitations for which transi- 
tion rates are available are subject to significant energy 
shifts caused by the quasiparticle interaction. For ijjijj 
and ipip* scattering states, these are effects of 0(iV -1 ) as 
discussed in Ref. |[ The exercise of caution is also indi- 
cated when the scaling behavior of the finite-size transi- 
tion rates changes at spectral boundaries as is frequently 
the case. 

Notwithstanding these caveats, the product ansatz is 
a useful tool for merging the best available transition- 
rate data and density-of-states data. It was successfully 
tested for the two-spinon excitations at M z = O.Eil Any 
major distortions of the lineshapes predicted by the prod- 
uct ansatz can be avoided if we omit all data points of 
M(q, ui) that are shifted across spectral boundaries. Any 
theoretical and computational advances that make it pos- 
sible to calculate transition rates for larger systems will 
improve the predictive power of the product ansatz. 



relation 



A. S D D(q,uj) and S zz (q,u>) 



S + -(q) = 2M z /N+S-+(q) 



(11) 



dictate s that the singularity is the one already described 



in Sec. 1VA : S+-(q) 



q 



l/tj-i 



The relative iptp* contribution to the integrated inten- 
sity is largest near the zone center and near the zone 
boundary as shown in the inset to Fig. 0(b). It gradually 
drops from 100% at q = to ~ 20% at the soft-mode 
position q — q s and then rises back to ~ 72.5% at q = it. 
Note that the TV-dependence of the relative intensity data 
is much stronger at q > q s than at q < q s . We shall see 
that the qualitatively different TV-dependences are also 
observed in transition rates, from which interesting con- 
clusions can be drawn. 



V. LINESHAPES 

To calculate the lineshapes of the ipip and i/ji/j* contri- 
butions to the dynamic spin and dimer structure factors 
we use, wherever applicable, the product ansatz 



S(q,w) = D(q,oj)M(q,uj) 



(12) 



discussed at some length in Ref. |fj. The factor D(q : uj) is 
the density of ipip or i^ip* states, which can be evaluated 



What can be observed in a fixed-g scan of the dimer 
and parallel spin fluctuations at q — q s l The lineshape 
determination of S^ D (yr/2,u)) is illustrated in Fig. |8| in 
comparison with corresponding data for Sf^ (tt/2, uj) as 
shown in Fig. 8 of Ref. ||. The first factor in the prod- 
uct ansatz is the same for both sets of data, namely the 
density ofipip* states in continuum P2. It has the charac- 
teristic shape with a square-root divergence at the upper 
band edge ljjj as shown in panel (a). 

The scaled transition rates Mj^ (tt/2,ui) and 
M^Jjr, (-71-/2, w) [panel (b)] are monotonically decreasing 
functions. For both kinds of fluctuations the data at low 
frequencies are consistent with the power-law divergence, 



Stf(*/2,u) 



(13) 



with 7/ — 2 = —0.468 ... as predicted by conformal in- 
variance. It is near ljjj where the two sets of data differ 
most. While both transition rate functions tend to van- 
ish at u>u, this tendency is considerably slower for the 
parallel spin fluctuations than for the dimer fluctuations. 
The resulting lineshapes are dramatically different. 

We saw that the slow approach to zero at ljjj of the par- 
allel spin transition rates Mf^ (tt/2,u)) combined with 
the divergence in (tt/2,ui) produces a diverging 
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FIG. 8: (a) Density of i/ji/j* states at q = ty/2 for N = 2048. 
(b) Dimer {DD) and parallel spin (zz) transition rates be- 
tween the psinon vacuum and the iptp* states at q — n/2 for 
N < 32 (zz) and JV < 28 (DD). (c) Lineshape at q = tt/2 of 
the t/>?/>* contribution to 5*22(9,0;) and Sooiq, w). All results 
are for M z = iV/4. The scales in (b) and (c) on the left (right) 
are for DD {zz). 



trend at lo v in Sff (7r/2,uA The result is a character- 
istic double-peak structure. El The more rapid approach 
to zero of the dimer transition rates M^ D (ir/2,ui) over- 
comes the divergence of (tt/2,oj) and produces a 
single-peak lineshape in Slp D (ir/2,u>). 

Now we consider the wave numbers halfway between 
the soft mode q s and the zone boundary or the zone cen- 
ter. The lineshape determination via product ansatz of 
S%% (q, w) at q = tt/4, 3tt/4 is illustrated in Fig. g. The 
corresponding data for the parallel spin fluctuations were 
shown in Fig. 9 of Ref. g. For this situation the upper 
edge of one band (g = tt/4) coincides with the lower edge 
of the other band (q = 37r/4). 

Even though the number of available data points for 
dimer transition rates is limited, there is a clear indica- 
tion that the lineshapes of the dimer fluctuations again 
consist of single-peak structures with a divergence at the 
lower band edge u>l 7^ and a shoulder reaching to the 
upper band edge Wp, in strong contrast to the double- 
peak structures predicted for the lineshapes of the paral- 
lel spin fluctuations. 

In summary, the spectral- weight distribution of the dy- 
namic structure factors which probe the dimer and par- 
allel spin fluctuations have many commonalities but also 
some very distinct properties. In both cases, the domi- 
nant spectrum is the continuum P2 of ijjtp* states within 
the invariant M z subspace which also contains the ground 
state \G). Both functions Slp D (q,w) and Sff {q,w) are 
strongly peaked along the lower continuum boundary 
(jJl{q)- Only the latter is also peaked along the upper 
continuum boundary tou{q). The divergence along wl(<z) 
is caused by the transition rates, whereas the divergence 
along uou{q) is a density-of-states effect. 



FIG. 9: (a) Density of inp* states at q = tt/4, 37t/4 for 
iV = 2048. (b) Dimer transition rates between the psinon 
vacuum and the tptf)* states at q — n/4,Sn/4 for TV = 16, 24. 
(c) Lineshape at q — 7r/4, 37r/4 of the VV'* contribution to 
S D D{q,u)- All results are for M z = N/4. 



B. S- + {q,cu) 

Here we focus on the lineshape at q = it of the ipi/j con- 
tribution to the dynamic spin structure factor S \.(q,u). 

The continuum P3 of tpij) states was previously found to 
be dominant. The results predicted on the basis of the 
product ansatz are shown in Fig. 

The density of states has the same characteristic shape 
as seen previously. The spectral weight is strongly 
concentrated in the lowest finite-size excitation. The 
scaling behavior of the transition rate for that state, 
A^|<G'|5'-|A)| 2 - (l/AO 1 /''- 2 , translates, via conformal 
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FIG. 10: (a) Density of ipip states at q = tt for N = 2048. 

(b) Perpendicular spin transition rates between the psinon 
vacuum and the ipip states at q = n for N = 12, 16, 20, 24, 28. 

(c) Lineshape at q = n of the tptjj contribution to S y{q,uj). 

All results are for M z = N/4. 
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invariance, into a power-law infrared divergence, 
Stt(n,u) ^uj 1 ^- 2 , I/77- 2 = -1.346... , 



(14) 



for the spectral-weight distribution of the infinite chain. 
The solid line represents a two-parameter fit, au; 1 ^ 1 ^ 2 + 
b, of the data points representing the lowest excitation 
for N = 12, 16, 20, 24, 28. The transition rate data at 
higher frequencies appear to approach zero sufficiently 
rapidly to overcome the divergent trend of the density 
of states to produce a monotonically decreasing spectral- 
weight distribution with a cusp singularity at the upper 
continuum boundary. 

Similar single-peak lineshapes are expected for fixed-g 
scans across the range of the continuum P3. Hence the 
ipip contribution to the perpendicular spin fluctuations 
is a structure that is strongly peaked along the lower 
continuum. boundary u)i,{q) in the shape of the psinon 
dispersiontj and a shoulder reaching to the upper bound- 
ary uJu{q) of the continuum P2. Given the strong di- 
vergence of S_A_(q,w) at ojL(q), the perpendicular spin 
fluctuations offer the most promising way to measure the 
energy-momentum relation of the psinon quasiparticle by 
means of neutron scattering. 



C. S+-(q,u) 

Here we are back to focusi ng on lineshapes produced 
by iptp* excitations as in Sec. V A , but not in the same 
invariant M z subspace. Nevertheless, the continuum P6 
as depicted in Fig. ^(a) produces, at q = tt/2, a band 
of equal width and location as continuum P2 depicted in 
Fig. I did. 

The data used in the product ansatz applied to 
S+_ (w/2,uj) are shown in Fig. ll[ The density of states 
is exactly the same as in Fig. |. The data for the tran- 
sition rates are monotonically increasing with lo. The 
trend in the low-frequency limit is that Afj__ (tt/2, lo) ap- 
proaches a finite value, possibly zero. Given the fact that 
jyi>-<i> u) is flat at low frequencies, the same trend is 

observed in Stz. (tt/2,u)). The prediction of conformal 
invariance for the leading infrared singularity is 



Sit (tt/2,lo)~^ 



2 = 3.531. 



(15) 



At the upper band edge lojj, the transition rate data 
exhibit a pronounced maximum which could either signal 
a divergence or a cusp singularity for the infinite system. 
Whatever the case may be, this enhancement amplifies 
the divergent density of states in Sr^L (n/2,ui). 

Recall that the parallel spin fluctuations exhibit a 
double-peak lincshape at the soft-mode wave number q s . 
By contrast, the lineshape of the perpendicular spin fluc- 
tuations at the soft-mode wave number q s is a single-peak 
structure with the spectral weight concentrated near the 
upper band edge. For other wave numbers we do at 
present not have enough transitions rate data points for 
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FIG. 11: (a) Density of ifiip states at q = n/2 for N = 
2048. (b) Perpendicular spin transition rates between the 
psinon vacuum and the ipip* states at q — n/2 for N = 
12, 16, 20, 24, 28. (c) Lineshape at q = tt/2 of the ipip* contri- 
bution to S+~(q,uj). All results are for M z = N/4. 



a useful application of the product ansatz. Nevertheless, 
from the few data points which we do have some inter- 
esting conclusions can be drawn. 

In the transition rate data currently available we ob- 
serve that the spectral weight is heavily concentrated in 
a single excitation for any given q 7^ 7r/2. For < q < q s 
that excitation is located along the lower boundary of 
continuum P6 and for q s < q < tx along the upper bound- 
ary. BothJpoundarics have the shape of the antipsinon 
dispersionJa Therefore, the contribution of Sfz. (g, u>) to 
the perpendicular spin fluctuations offers the most direct 
way to measure the dispersion of the antipsinon quasi- 
particle by means of neutron scattering - not once but 
twice, in different parts of the Brillouin zone. 

The outstanding role of the excitations along the two 
antipsinon branches at q < q s and q > q s is illustrated 
in Fig. [12], where we have plotted the transition rates 
KGIS+IA)! 2 versus q of all states from continuum P6 
across various system sizes. However, the role of the 
finite-TV excitations which are part of these dominant 
branches is different at q < q s and q > q s . This is ev- 
ident by comparison of panels (a) and (b) which show 
differently scaled transition rates N a \(G\S+\\)\ 2 . 

To make sense of the data in this representation we 
must distinguish three scenarios and ask the question: 
For what value of the scaling exponent a do the transition 
rate data exhibit minimal TV-dependence? 

(i) For states that are inside a continuum, the product 
ansatz requires requires that the exponent a = 1 mini- 
mizes the A^-dependence of the transition rate data. 

(ii) For states that mark the boundary of a continuum 
where the spectral-weight distribution has a divergent 
singularity that is not caused by a divergent density of 
states the exponent which minimizes the A^-dependence 
of the transition rate data is in the range < a < 1. 
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FIG. 12: (a) Unsealed and (b) scaled iptfi* transition rates for 
S+- (q, to) at M z = N/4 for AT = 12, 16, 20, 24, 28. 



(iii) For states that form a branch with nonvanishing 
spectral weight in the limit N — > oo, the exponent that 
minimizes the TV-dependence is a = 0. 

There is ample evidence for the first scenario in the 
results presented earlier. The data in panel (a) strongly 
indicate that the third scenario is realized for the an- 
tipsinon branch at < q < q s . The evidence is rigorous 
for the isolated excitation at q — which carries all the 
spectral of (q,w) as discussed previously. 

Panel (b) suggests that the second scenario applies to 
the antipsinon branch at q s < q < tt. The exponent 
a = l/rj = 0.653. . . used here is suggested by confor- 
mal invariance, but there is a strong possibility that the 
singularity of S+_ (q, lo) along the antipsinon branch at 
Qs < <Z < t is gov erned by a g-dependent exponent. The 
insets to Fig. [12] show the iV-dependence of the corre- 
sponding data at q = tt, which is off the scale used in the 
main plots. At the zone boundary, the width of contin- 
uum P6 has shrunk to zero, which is likely to affect the 
./V-dependence of the transition rate singularly, as indi- 
cated. 

In conclusion, this study makes predictions of unprece- 
dented detail for the lineshapes of the spin and dimcr 
fluctuations in the ID s = -z Heisenberg antiferromagnet 
at zero temperature and nonzero magnetic field. Fore- 
most among these predictions is the direct observability 
in the perpendicular spin fluctuations of the dispersion 
relations for the two quasiparticles which play a crucial 
role in this situation: the psinon and the antipsinon. 
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APPENDIX A: CALCULATING MATRIX 
ELEMENTS VIA BETHE ANSATZ 

The Bethe ansatzi is an exact method for the calcu- 
lation of eigenvectors of integrable quantum many-body 
systems. The Bethe wave function of any eigenstate of 
(|l|) in the invariant subspace with r = N/2 — M z reversed 
spins relative to the magnon vacuum, 

W= E a(n 1 ,...,n r )S- i ---S- r \F), (Al) 

l<ni<...<n r <A r 

has coefficients of the form 



o(m, . .. ,n r ) = exp li^k-pjTij + - ^6 VtVj 



(A2) 



determined by r magnon momenta fcj and one phase angle 
9ij = —Oji for each magnon pair. The sum V G S r is over 
the permutations of the labels {1,2,... , r}. 

The consistency requirements for the coefficients 
a(ni,... ,n r ) inferred from the eigenvalue equation 
H\ip) — E\ip) and the requirements imposed by transla- 
tional invariance lead to a set of coupled nonlinear equa- 
tions for the ki and 6ij. A computationally convenient 
rendition of the Bethe ansatz equations for a state spec- 
ified by Bethe quantum numbers I\ , . . . , I r has the form 



N<p(zi) = 2-Kh 



Zj )/2] 



1, 



(A3) 



where 4>(z) = 2arctanz, ki — tt ~ 4>{zi) and 6ij — 
7rsgn[3?(zj — Zj)] — 4>[(zi — Zj)/2] . The energy and wave 
number of the eigenvector thus determined are 



E-E f 



^ 1 

i=i 



TXT 



N ^ 1 

i=l 



(A4) 



where Ep — JN /A is the energy of the magnon vacuum. 

In the past, The Bethe ansatz was rarely used for the 
purpose of calculating matrix elements. The main de- 
terrent has been the need of evaluating the sum V 6 S r 
over the r! magnon permutations in the coefficients ( |A2| ) 
of the Bethe eigenvectors (Al). However, the tide is now 
changing rapidly for two reasons: (i) the availability of 
vastly higher computational power, (ii) theoretical ad- 
vances that make it possible to reduce matrix elements 
of Bethe wave functions to determinantal expressions $3 

In the following, we sketch how the matrix elements 
can be manipulated effectively by using the Bethe wave 
functions directly. The use of the determinantal expres- 
sions of Bethe ansatz transition rates for the calculation, 
of dynamic structure factors will be reported elsewhere.EEl 

In designing an efficient algorithm, we must heed the 
fact that in the calculation of a single matrix element, 
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the sum V £ S r is evaluated many times, once for ev- 
ery coefficient a(n\, ■ ■ ■ ,n r ) of the two eigenvectors in- 
volved. Under these circumstances, it is imperative that 
the algorithm has rapid access to a table of permuta- 
tions. Such tables, can be generated recursively by pow- 
erful algorithms.Eil 

The computational effort can be redu ced considerably 
if we use the translational symmetry of (Al), guaranteed 



by the relation a{n\ + . 







a(n±, ... ,n r ) 



between Bethe coefficients (A2) pertaining to basis vec- 
tors that transform into each other under translation. 
Translationally invariant basis vectors have the form 



\r,k) 



i 



E 



(A5) 



•o 1=0 



where the a,j = a{n^\ . . . ,ni- J '), the Bethe coefficients 
of the generator basis vectors |j)o, are the onl y on es that 
must be evaluated. The last expression of (A7) holds 



,0'h 



because the Bethe coefficients a,j of all generators \j)o 
which do not occur in the set are zero. 

We calculate transition rates for the dynamic structure 
factor (|^) in the form 



\(G\SZ\m)f 



Uo\\ 2 \\4> m \\ 2 



(A8) 



where \ifio), \ip m ) are the (non-normalized) Bethe eigen- 
vectors of the ground state and of one of the excited states 
from classes (i)-(vi), respectively. The norms are 



where \j)i = T l \j) Q 



I. 



, n 



1} and 1 < 



N/dj < N is an integer. The wave numbers k realized in 



the set (A5) are multiples mod(27r) of 2-K/dj. 

The set of basis vectors \j) = , ... , rir ),j 
1, . . . ,d, are the generators of the translationally invari- 
ant basis. The set of distinct vectors \j; k) for fixed k 
is labeled j € Jk C {1, . . . , d}. The rotationally invari- 
ant subspace for fixed N/2 — r, which has dimensionality 
D = N\/[r\(N — r)!], splits into N translationally invari- 
ant subspaces of dimensionality one for each wave 
number k — 2m/N, n = 0, . . . , N — 1. We have 



= 5>Kf- (A9) 
j'=i 

The matrix element (ipo\S£\ip m ) is nonzero only if q = 
k m — ko + 27rZ. For the fluctuation operator S* it can be 
evaluated in the form: 



{tpo\Sq\lp m )- 



d N 

1 Y^(°)„M' 



dt-X 



3 3 t[ to 



j lq i(j\S z n \j)i. 
(A10) 



D = Y, d i= E D ^ D *= E 

j=i o<fc<27r jeJk 



(A6) 



The Bethe eigenvector (Al) expanded in this basis can 
thus be written in the form 

dj — l d dj—x 

iv^> = E °i E ellk \^ = E a i E ellk \^ ( A7 ) 

jeJk i=o j=x i=o 



The non- vanishing matrix elements (ipo\S^\ipm} for the 
other spin fluctuation operators must also satisfy q — 
k m — ko+2'K t Z and can be reduced to somewhat more com- 
plicated expressions involving elements i (jo\S^\j m )i m 
between basis vectors from different Sj. subspaces. 

The memory requirements for the calculation of one 
such matrix element are 6.7MB for N — 18, r = 9 and 
73MB for N = 20, r = 10. 
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